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We apply recently developed 1/{N — 1) expansion to a particle-hole asymmetric SU(A^) Anderson 
model with finite Coulomb interaction U. To leading order in l/( A'^ — 1) it describes the Hartree-Fock 
random phase approximation (HF-RPA), and the higher-order corrections describe systematically 
the fluctuations beyond the HF-RPA. We show that the next-leading order results of the renormal- 
ized parameters for the local Fermi-liquid state agree closely with the numerical renormalization 
group results at A = 4. It ensures the reliability of the next-leading order results for A*' > 4, and 
we examine the N dependence of the local Fermi-liquid parameters. Our expansion scheme uses the 
standard Feynman diagrams, and has wide potential applications. 
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I. INTRODUCTION 

Exploring the ground state and the low-energy exci- 
tations of strongly correlated electron systems is one of 
the fundamental subjects of quantum many-body theory. 
The low-lying energy states determine not only equilib- 
rium properties but nonequilibrium transport deviating 
from the linear response. For instance, in quantum dots 
the universal Kondo behavior of the steady currentii"— 
and shot noise^"— at small bias voltages can be described 
in terms of the quasi-particles of a local Fermi liquid. 

The universal behavior has been studied in more gen- 
eral situations with the orbital degeneracy.—"— The es- 
sential feature of the universality in this case may be 
deduced from the low-energy states of the A^-fold degen- 
erate Anderson impurity. The exact numerical renormal- 
ization group (NRG) approacbi^ was successfully applied 
to this model for small degeneracies N < A^^"— which 
for N = 2 corresponds to the spin degeneracy. However, 
alternative approaches arc required for large degeneracies 

> 4. 

The conventional 1/N expansion, or non-crossing ap- 
proximation (NCA), is one of such approaches 1^^*^^ This 
and related approximations are based on the perturba- 
tion theory with respect to the hybridization matrix el- 
ement V between the impurity and conduction band. 
They treat mainly the electron filling less than one elec- 
tron in A^ orbitals, assuming an infinite Coulomb inter- 
action U. Extensions of the methods of this category 
to the case of finite U have been developed, also on the 
basis of the perturbation theory in V, using the slave- 
boson or the resolvent technique.— However, in order 
to explore a wide range of the electron fillings and of the 
other parameters, it is still necessary to develop different 
approaches applicable to the low-energy excitations. 

In a previous paper, we proposed a completely different 
approach, using a scaling that takes u = (A^ ^ as an 
independent variable instead of the Coulomb interaction 
U alone.— The factor A^ — 1 appears naturally as it rep- 
resents the number of interacting orbitals, excluding the 
one prohibited by the Pauli principle. With this scaling 



the perturbation series in U can be reorganized as an ex- 
pansion in powers of 1/{N — 1). This can be carried out 
by using a standard Feynman diagrammatic technique 
for the fermions with two-body interactions. Specifically, 
diagrammatic classification, similar to the one for the 
A^-component cp'^ model for critical phenomena,— is ap- 
plicable with an extension that takes into account the 
fermionic anticommutation relation. We calculated the 
renormalization parameters for the quasi-particles up to 
order 1/(A^ — 1)^ at half-filling, and confirmed that the 
results agree closely with the NRG results for A' = 4. It 
enabled us to clarify the A^ dependence of the nonequi- 
librium Kondo behavior for A^ > 4;^^ 

In the present paper, we apply this approach to away 
from half-filling, and examine the ground-state proper- 
ties. To this end, we choose the unperturbed Green's 
function such that it includes the level shift due to the 
a; = component of the self-energy, which at zeroth order 
is given by the Hartree-Fock (HF) approximation. The 
first order contributions in the 1/(A^ — 1) expansion con- 
sists of the vertex and self-energy corrections, which are 
equivalent to those of the random phase approximation 
(RPA);^ The higher order contributions systematically 
capture the correlation effects beyond the HF-RPA. We 
show that the order 1/(A^ — 1)^ results of the renormal- 
ized parameters agree well to those of the NRG also in 
the particle-hole asymmetric case. It ensures an early 
convergence of the 1/(A^ — 1) expansion for A^ > 4. With 
this approach, we also examine the A" dependence of the 
local Fermi-liquid parameters. Our expansion procedure, 
using the standard Feynman diagrams, is quite general. 
It can be applied to the Keldysh formalism, and to lattice 
systems such as the Hubbard model. 

This paper is organized as follows. We describe the 
model and the basic idea of the scaling of the Coulomb 
interaction U in Sec. [Hi The diagrammatic classification 
for the leading and next-leading order terms in 1/{N— 1) 
expansion is explained in Sec. IIIII The numerical results 
for the renormalized parameters are presented in Sec. lIVI 
A brief summary is given in Sec. |Vl 
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II. MODEL AND FORMULATION 

We consider the SU(Af) impurity Anderson model with 
finite interaction U , given by "H = Hd+'Hc + 'H'V + 'Hu, 
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Here, d}^ is the creation operator for an electron with 
energy and orbital m (= 1,2,--- , A^) in the impu- 
rity site, and n^^^^ = rfj„d„j. The operator c|„j creates a 
conduction electron that is normalized as {C(:m,c^,^,} = 
5{e— e'). The hybridization energy scale is given by 
TTpV^ with p = 1/{2D). We consider the parameter 
region where max(A, |erf|, U) D. 

We use the imaginary-frequency Green's function that 
takes the form for \uj\ <C I?, 
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where S](ia;) is the self-energy due to T-lu. The ground- 
state average of the local charge can be deduced from 
the phase shift 5 = cot-'^{E*/A) with E* = Ed + E(0), 
using the Friedel sum rule (ndm) = S/n. The deriva- 
tive of T,{iuj) determines the renormalized parameters 
7=1- d'S{iuj)/d{iuj)\^^o, z = 1/7, = ^ E^, and 
A = zA. Furthermore, using the Ward identities;^ the 
enhancement factor for spin and charge susceptibilities, 

Xs = Xmm ~ Xmm' Xc = Xmm + ( ~ 1 ) Xmm/ ' '"^^ 

expressed in terms of the renormalization factor 7 and the 
vertex function r^,„,.,„,^(iiui, ia;2; iws, ioj^) for m ^ to'. 
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The residual interaction between the quasi-particles is 
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We showed in the previous work that it is essential 
for > 2 to scale the Coulomb interaction by multiply- 
ing a factor A' — 1 in the particle-hole symmetric casejii 
Away from half-filling, it is more natural to include an 
additional factor sin^S such that 
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as the local density of states at the impurity site is given 
by -ilmG(iO+) = sin2,5/(7rA). Similarly, we introduce 
the renormalized version of the scaled coupling constant. 
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The Wilson ratio R and that for the charge sector can be 
expressed in terms of the renormalized coupling AT, 
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Conversely, it can also be written as AT = (A' — l)(i? — 1). 
By definition, K is bounded in the range < AT < 1, 
and approaches A" — >■ 1 when the charge fluctuation is 
suppressed Xc 0. At half-fiUing ed = -(A^-l)C//2 and 
Ed = 0, it has been confirmed that K varies in a narrow 
region as A^ increases, and converges rapidly to the RPA 
value that is asymptotically exact in the A^ — > 00 limit.— 
Away from half-filling, the renormalized impurity level 
E^ determines the charge distribution near the impurity, 
and should be calculated in an optimal way. We treat 
this problem by taking the unperturbed Green's function 
such that 
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In this set up, the energy shift due to S(0) is included 
in the unperturbed part of the Hamiltonian, defined by 
'Ho = 'Hd + 'Hc + Hy + X]m^(0) "^dm , ^^nd it is subtracted 
in the perturbation Hamiltonian 
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Here, the counter term is defined formally by A = E(0).— 
Then, the full Green's function defined in Ec^. ^ can be 
rewritten in the form 
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The redefined self-energy Yi(iui) represents the correc- 
tions due to T-Lui and can be calculated perturbatively 
as a function of E^ and A. These two parameters can be 
determined through the renormalization condition 
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and Eq. (j4]) that can be rewritten, using the Friedel sum 
rule, as 
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where ^d = ed + (N - l)U/2. At half-filhng = 0, both 
of the two parameters vanish, A^ = and A = 0. 



III. 1/(A - 1) EXPANSION 

The perturbation expansion with respect to can 
be classified according to the power of 1/(A^ — 1), by 



FIG. 1. The leading order diagrams in the 1/{N — 1) ex- 
pansion. The wavy and sohd hnes indicate the Coulomb re- 
pulsion U and unperturbed Green's function Go, respectively. 
The double wavy line represents the sum of the bubble dia- 
grams, which corresponds to Uhuhii^) given in Eq. (|17|) . Each 
fermion loop gives a factor of order (A'^ — 1) through a sum 
over the different orbitals [See Appendix \K\ . 



choosing the scaled interaction g as an independent vari- 
able instead of bare U . The counter term can also be 
expanded such that 




FIG. 2. The diagrams which provide the order 1/{N — f) 
contributions with higher order corrections. The cross denotes 
— A of the counter term. 



Fig. [21 which determines the leading order term of K 
with the HF value of K* , as 
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and the coefficient Afc is determined by the requirement 
of the renormalization condition Eq. (|13p is satisfied to 
each order in the 1/{N — 1) expansion. 



A. Zeroth order 

By construction, the tadpole diagrams and the contri- 
butions arising from the bilinear part of Hu defined in 
Eq. ©, i.e., -UY^^_^^,{ndm') ndrn, cancel each other. 
Thus at zero order, the counter term vanishes A = 0, and 
is determined by Eq. (fT4|) . which in this case gives 
the HF value, 
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Note that the corrections due to the RPA are still not 
included in the zeroth order. 



B. Leading order 

The leading order corrections in the 1/{N — 1) expan- 
sion arise form a series of the bubble diagrams of the 
RPA type indicated in Fig. [U and can be expressed in 
the form 
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The precise form of lA-^^-^{iu) and XQ(iw) are provided 
in Appendix |X1 The propagator U^^'^^{iu!) contains also 
the higher order contributions in the ^/{N — 1) expan- 
sion. This is because the orbital indices for adjacent bub- 
bles have to be different due to the Pauli exclusion, and 
thus the summations over internal m's for the fermion 
loops arc not independent. The order 1/{N — 1) terms 



of r. 
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(0, 0; 0, 0) arises from the first diagram in 



Note that 7rAxo(0) = 1/[1 + {E*/A)'^] also equals sin^S. 
Equation ((T9| determines the Wilson ratio R, defined in 
Eq. dHI, to order 1/{N - 1). _ 

The leading order self-energy due toJiu arises from the 
second and third diagrams in Fig. [2] These two diagrams 
represent the contributions of the RPA fiuctuation and 
counter term, and can be expressed in the form 
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where A^' = — 1. Therefore, from the renormalization 
condition defined in Eq. (jl3p the counter term can be 
deduced to the first order in 1/(A^ — 1), 

A = S(V^')(0) + 0(^^) . (22) 

In a similar way the coefficient Afc of the counter term 
can be determined order by order. We can also calculate 
the leading order term of 7 from E*^^/^ 
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Here, j^lP^' denotes the next-leading order contribu- 
tions, which will be taken into account later. 



C. Next leading order 

Fluctuations beyond the RPA enter through the next- 
leading and higher order terms in the 1/(A^ ~ 1) ex- 
pansion. In the particle-hole asymmetric case, the or- 
der 1/{N— 1)2 contributions to r^,^,.^,^(0, 0; 0, 0) arise 
from the diagrams shown in Figs. [3] and HI and also from 
the order 1/{N— 1)^ order component oi U^^^^^{iu!) shown 
in the first diagram of Fig. ^ i.e., the second term in 




FIG. 3. The order l/{N-lf dia fframs for the vertex function 
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(0,0; 0,0) for m / m'. 







FIG. 4. The order 1/{N — 1)^ vertex diagrams for m ^ m' 
contribute away from half-fiUing in addition to the ones given 
in Fig. [3l Tire contribution of the particle-particle pairs and 
that of particle-hole pairs cancel each other at half-filling. 




FIG. 5. The order l/{N-lf self-energy diagrams contribute 
to z (= 1/7). 




FIG. 6. The order l/{N-lf self-energy diagrams contribute 
away from half-filling in addition to the ones given in Fig. 
[5] The contribution of the particle-particle pairs and that 
of particle-hole pairs cancel each other in the particle-hole 
symmetric case. The counter term also vanishes, A = 0, at 
half-filling. 



the right-hand side of Eq. ([T7)) . Note that the contri- 
butions of the diagrams shown in Fig. |4] vanish at half- 
filling. This is because the particle-particle and particle- 
hole pairs give opposite contributions, and cancel each 
other in the particle-hole symmetric case. 

Summing up the contributions from all these diagrams, 
we calculate the Wilson ratio i?, defined in Eq. ([3]), to 
order l/(iV— 1)^. The result of the renormalized coupling 
K can be expressed in the following form, which is exact 
up to terms of order 1/{N — 1), 
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In this formula the contribution described by I^, emerg- 
ing in the numerator in the right-hand side of Eq. (|25p , 
becomes finite in the particle-hole asymmetric case, and 
it represents the contributions from the eight diagrams in 
Fig. 21 excluding the first two caused by the counter term. 
The prefactor K* /{I + K*) in the right-hand side of Eq. 
(p5)) is required to be calculated to order l/(iV — 1), by 
including the energy shift due to Eq. (HI]). Alternatively, 
through Eq. pi]) with Eq. ([H]), the bare parameter 
can be determined as a function of E*^ to order 1/{N— 1). 

We also consider the next-leading order term of the 
self-energy, which we found in the previous work gives 
important contributions to the renormalization factor z 
in the particle- hole symmetric caseJ^ The order 1/{N — 



1)^ self-energy corrections arise from the diagrams shown 
in Figs. [5] and [6l Furthermore, the corrections arise also 
from the second and third diagrams in Fig. [2] through 
the higher-order component of Z//bub(iw) and that of A, 
respectively. In the particle-hole asymmetric case, the 
diagrams indicated in Fig. |6] give finite contributions as 
the exact cancellation between the particle-particle and 
particle-hole pairs occurs only at half-filling. We have 
calculated all these next-leading order contributions of 
the self-energy numerically to obtain 7, A, and the 
phase shift i5 to order 1/(7V — 1)^. 



IV. NUMERICAL RESULTS 

We have also carried out the NRG calculations for 
= 4, and have compared with the results of the 
1/(A — 1) expansion. We find generally that the next- 
leading order results agree well with the NRG results, 
especially for small interactions g < 3.0. Typical exam- 
ples are presented in Figs. [7] and [H] 

In Fig. [71 comparisons are made for {ridm), sin^S and 
Erf, choosing a rather large value g = 9.0 for the interac- 
tion in order to show clearly the small deviations. We see 
the very close agreement between the order 1/(A^ — 1)^ 
results (red sold line) and the NRG (•) results over the 
whole region of £,d- The order 1/(A^ — 1) results, cor- 
responding to the HF-RPA, are also plotted in Fig. [7] 
with the green dash-dot line although they are almost 
concealed under the red solid line. Nevertheless, a very 
small deviation is visible near ^^(/(TrAg) ~ ±0.4, showing 
that the order 1/(A^— 1)^ corrections improve the results 
slightly closer to the exact NRG results. It indicates that 
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FIG. 7. (Color online) (ndm) ~ S/tt and s'm^S (upper panel), 
and (lower panel) plotted as a function of for A'^ = 4 
and g = 9.0. The circles (•) represent the NRG results, and 
the red solid line the order 1/{N — 1)^ results. The blue 
dashed line denotes the zeroth order results obtained with the 
HF approximation without the RPA corrections. The order 
1/{N — 1) results, corresponding to the HF-RPA, are also 
plotted with the green dash-dot line although it is almost 
concealed under the red solid line. 



the phase shift S can be approximated reasonably by the 
order l/(iV— 1) self-energy given in Eq. ([?!]) . which signif- 
icantly improves the simple HF values (blue dashed line) 
given at the zeroth order without taking into account the 
RPA fluctuations. We have also confirmed similar trend 
for a larger value of the repulsion g = 12.0. We also see 
in the lower panel of Fig. [7] that the rcnormalizcd impu- 
rity level stays close to the Fermi level in the region 
|^d|/(7rA(7) < 0.3 as the interaction g is relatively large 
in this case. Outside of this region the impurity states 
tend to be empty at ^^/(TrAg) > 0.5, or fully occupied 
at £,d/{TrAg) < -0.5. We note that the NRG resuhs 
for sin^ 5 and {ndm) wave weakly at ^(^/(TrAg) ~ ±0.3, 
relating to the Coulomb oscillation discussed below. 

In Fig. [51 the renormalized parameters K and z are 
plotted for = 4, choosing g = 2.0 and 9.0. We see that 
the next-leading order corrections (red solid line) signifi- 
cantly improve the leading order results corresponding to 
the HF-RPA (blue dashed line), towards the exact NRG 
values (solid circles). Specifically, for small interactions 
{g = 2.0), the close agreement between the next-leading 
order results and that of the NRG are seen both for K 
and z. However, for large interactions {g = 9.0), the NRG 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

FIG. 8. (Color online) Renormalized parameters K (upper 
panel) and z (lower panel) plotted as a function of for 
= 4, choosing g — 2.0 and 9.0. The circles (•) represent 
the NRG results. The blue dashed line represents the leading 
order results, corresponding to the HF-RPA. The read solid 
line represents the next-leading order results, for which the 
fluctuation beyond the RPA are taken into account. 



results show an oscillatory behavior, as reported recently 
also in Ref. This behavior is seen more clearly for 
the renormalization factor z than the coupling K. The 
oscillation is caused by addition of electron into the lo- 
cal orbitals, which in the atomic limit of A = occurs 
discontinuously at —Cd/U = 0, 1, 2, . . . , — 1. Thus, for 
strong Coulomb interactions, iV — 1 local minima (max- 
ima) appear for z (K) in the region of [^^[/(TrAg) < 0.5, 
which corresponds to — (iV — 1)U < ed < 0. 

Truncation of the 1/{N — 1) expansion does not re- 
produce the oscillatory behavior. Nevertheless, we see in 
Fig. [8] that the next-leading order results of K are still 
close to the exact NRG results even for large g = 9.0, 
especially at the plateau region for |^t;|/(7rA(7) < 0.2, 
where K reaches near the unitary-limit value 1.0 and 
the charge excitation is suppressed significantly. Further- 
more, we also see that the order 1/{N — 1)^ results of z 
for g = 9.0 approach closely the maxima, emerging at 
^((/(ttA^) ~ ±0.17 in the NRG results, which correspond 
to mixed-valence regions in between two adjacent inte- 
ger charge states. Therefore, the fluctuations beyond the 
RPA, taken into account to order 1/{N — 1)^, can de- 
scribe approximately an upper envelop for the waving 
curve of the renormalization factor, z, for A^ > 4. Note 
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FIG. 9. (Color online) Next leading-order results plotted 
vs for g = 9.0, for several {— 4,6, . . . , oo): K and z are 
calculated up to order 1 / (A'' — 1) and 1 / (A'^ — 1)^ , respectively. 
The orange dash-dot line represents the NRG results for N = 
2, and the gray dashed line represents the exact results in the 
N ^ oo limit. 
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FIG. 10. (Color online) Next leading-order results plotted 
vs for g = 2.0, for several (= 4, 6, . . . , oo): K and z are 
calculated up to order 1 / (A' — 1) and 1 / (A'^ — 1)^ , respectively. 
The orange dash-dot line represents the NRG results for A' = 
2, and the gray dashed line represents the exact results in the 
TV — i- oo limit. 



that the Kondo energy scale can also be estimated from 
this result^ of z as the width of the Kondo resonance is 
given by A = zA. 

In order to clarify the dependence of the renormal- 
ized parameters on the orbital degeneracy N, the next- 
leading order results of K and z for several values of N 
(= 4, 6, . . .) are shown in Figs. 1^ and [TUl for g = 9.0 and 
2.0. In these figures, the NRG results for = 2 (or- 
ange dash-dot line) and the exact values in the — > oo 
limit (gray dashed line) are also shown. We see, in both 
strong g = 9.0 and weak ^ = 2.0 coupling cases, the value 
that the renormalized coupling K can take at > 4 is 
bounded in a narrow region between the curve for = 4 
and the one for the large N limit that is described by 
the RPA given in Eq. ([TO]). As N increases, K converges 
rapidly to the asymptotic value in the N ^ oo limit. 

The zeroth-order HF results for {ridm), sivl'S and e^, 
which are plotted in Fig. [7] with the blue dashed lines, 
also correspond to the exact values in the large N limit. 
Therefore, as N increases from = 4 to oo, the phase 
shift 5 and the renormalized impurity level vary also 
in a narrow region between the red-solid and blue-dashed 
lines indicated in Fig. [T] 

In contrast, we sec in the lower panel of Figs. [S] and 
1101 that the rcnormalization factor z varies in a wider 
range. Specifically, for large interaction g = 9.0, it varies 



in the range of 0.4 < 2 < 1.0 at |^d|/(7rA.g) < 0.3 as the 
orbitals increases from A^ = 4 to A^ — > oo. As mentioned 
for Fig. [51 the next-leading order results of z become less 
accurate than that of K for large g. Nevertheless, we 
also see in Fig. [8] that for small interactions {g = 2.0) 
the order 1/(A^ — 1)^ results of z are very close to the 
exact NRG values even for A^ = 4 where the degeneracy 
A^ is not so large. Furthermore, as A^ increases, the order 
1/(A^— 1)^ corrections dominate the fluctuations beyond 
the RPA, and capture reasonably the correlation effects, 
especially in the mixed-valence regime. 

V. SUMMARY 

In summary, we have described the 1/{N — 1) expan- 
sion approach to the SU(A^) Anderson model with finite 
Coulomb interaction, without assuming the particle-hole 
symmetry. We have calculated the renormalized param- 
eters up to order 1/(A^ — 1)^. Specifically, the results 
of the phase shift S, or E^, show very close agreement 
with the NRG results for A^ = 4. This trend has been 
confirmed for g < 10.0, from small to relatively large in- 
teractions. The scaled Wilson ratio K ~ {N — l)(i? — 1) 
and the rcnormalization factor z also agree well with the 
exact value in the mixed-valence regime, or not too large 
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interactions. As N increases from N = A, K varies in 
a narrow range and converges rapidly to the RPA value 
that is asymptotically exact in the large N limit, whereas 
z varies in a wider range. This behavior of K and z is de- 
termined essentially by the fluctuations beyond the HF- 
RPA, captured through the order 1/{N— 1)^ corrections. 

The 1/{N ~ 1) expansion provides a well-defined and 
controlled way to take into account the fluctuations be- 
yond the HF-RPA. Furthermore, the expansion scheme 
based on the standard Feynman diagrammatic formalism 
for fermions with two-body interactions is quite general 
and flexible. Therefore, this approach has wide potential 
applications, such as the nonequilibrium problem in the 
Keldysh formalism and the Hubbard-type lattice models. 
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Appendix A: Coefficients for the bubble diagrams 

The contribution of the bubble diagram on the vertex 
correction for m ^ m' , shown in Fig. [U is not a simple 



geometric series in powers of U, 
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where xoii^) represents the contribution of a single loop, 
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Specifically, 7rAxo(0) = 1/[1 + (-5^/^)^] in the static 
limit w — >■ 0. The coefficient Ak in Eq. (|Aip arises from 
the summations over the orbital indices for the series of 
k fermion loops, which can be calculated as 



mi m2 fUk 



(7V-l)fe+i-(-l)fe+i 



N 



p=0 



-1 
N -1 



(A3) 



Here, the primed sums for m^'s run under the restriction 
due to the Pauli exclusion; rrij ^ nij+i for j = 0, 1, . . . , fc 
with Too = TO and rrik+i = to'. Equation (jA3[) shows 
that the coefficient Ak can be regarded as a polynomial 
of (A^— 1), and Z-/bub(*w) has the higher-order components 
in the expansion with respect to 1/(A^ — 1). 
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